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SUMMARY 


An explicit time/space finite difference procedure is used to model the propaga- 
tion of sound in a quasi one-dimensional duct containing high Mach number subsonic 
flow. Nonlinear acoustic equations are derived by perturbing the time-dependent 
Euler equations about a steady, compressible mean flow. The governing difference 
relations are based on a fourth-order, two-step (predictor-corrector) MacCormack 
scheme. Difference equations for the source and termination boundary conditions are 
derived from the appropriate characteristic relationships. The solution algorithm 
functions by switching on a time harmonic source and allowing the difference equa- 
tions to iterate to a steady state. A significant advantage with this approach is 
that the nonlinear terms can be retained and evaluated with only modest additional 
computer cost above that required for a linear model calculation. 

The principal effect of the nonlinearities was to shift acoustical energy to 
higher harmonics. With increased source strengths, wave steepening was observed. 

This phenomenon suggests that the acoustical response may approach a shock behavior 
at a higher sound pressure level as the throat Mach number approaches unity. Where 
applicable, comparisons were made with the calculations from a linear finite element 
algorithm. On a peak level basis, good agreement between the nonlinear finite dif- 
ference and linear finite element solutions was observed, even though a peak sound 
pressure level of about 150 dB occurred in the throat region. Nonlinear steady state 
waveform solutions are shown to be in excellent agreement with a nonlinear asymptotic 
theory. 


INTRODUCTION 

It has been customary in computational duct acoustics to base calculative pro- 
cedures on steady state linear equations using direct implicit numerical schemes. 

The Galerkin/weighted-residual method summarized in reference 1 is a state-of-the-art 
example of this approach. Such direct methods, whether based on finite element or 
finite difference methodology (ref. 2), produce a large set of coupled linear equa- 
tions, which are usually solved by matrix manipulations, that is, inversion/ 
decomposition procedures. However, many practical aeroengine duct modeling situa- 
tions require extremely large numbers of degrees-of-f reedom for accurate resolution. 
Not only are the direct methods costly in computer time for these problems, but in 
many instances, the computer core storage limitations are exceeded, and expensive 
out-of-core solution algorithms are necessary. 

Recently, Baumeister (refs. 3 to 5) advocated the use of explicit time-dependent 
difference methods for steady state duct acoustic calculations. With these methods, 
the time-dependent equations, subject to a harmonically varying source boundary con- 
dition, are iterated over a time/space difference discretization grid until steady 
state is achieved. Such procedures have the advantage of tremendously reducing core 
storage requirements. In addition, both linear and nonlinear equations can be 
handled, since the nonlinear terms are simply retained in the difference relations 
and evaluated during each iterative calculation step with a minimum of additional 
cost. 



In this paper, an explicit time/space finite difference procedure is used to 
simulate the propagation of sound in a quasi one-dimensional variable area duct con- 
taining high Mach number subsonic flow. For this physical situation it has been 
shown, both analytically and experimentally (refs. 6 to 8, for example), that signif- 
icant intensification of the acoustic wave structure occurs in the throat region as 
the flow velocity approaches sonic value. In fact, linear theory is singular for a 
throat Mach number of unity (ref, 6) and is therefore invalid in the throat region. 
Nonlinear acoustic equations are derived by perturbing the time— dependent Euler equa- 
tions about a steady compressible state. The acoustic equations are then cast into a 
nondimensional form and discretized in space and time by a fourth-order, two-step 
(predictor-corrector) MacCormack difference scheme (ref. 9). Numerical boundary 
conditions are derived for the source and termination sections from the appropriate 
characteristic relationships. The solution algorithm functions by switching on a 
real time harmonic source and allowing the difference equations to iterate to a 
steady state condition. 

Calculated results are presented for several area reduction ratios, throat Mach 
numbers, and acoustic source sound pressure levels. Where applicable, comparisons 
are made with the results calculated by a linear Galerkin/weighted-residual finite 
element algorithm. At higher sound pressure levels and throat Mach numbers approach- 
ing unity, for which nonlinear effects become significant, results calculated by the 
finite difference procedure are compared with the asymptotic theory developed by 
Callegari and Myers (ref. 10). 


SYMBOLS 

A area variation 

c speed of sound 

f frequency 

J number of grid points 

j node index, x = (j - 1 ) Ax 

k wave number, k = 

L length of flow duct 

M Mach number 

N number of time steps 

n time index, t = n At 

p acoustic pressure 

S source amplitude 

T^,T 2 characteristic quantities, equations (19) and (21) 
t time 
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u acoustic particle velocity 

u acoustic streaming velocity 

X axial coordinate 

Y ratio of specific heats 

Ax = L/(J - 1) 

z perturbation factor 

p acoustic density 

(0 circular frequency, co = 2nf 

Subscripts: 

j quantity evaluated at jth node located at x = (j - 1) Ax 

o ambient property 

s mean flow property 

t throat quantity 

Superscripts : 

n quantity evaluated at time t = n At 

(1) predictor quantity 

A tilde (~) over a symbol indicates a dimensional variable, A caret ('") over a 
symbol indicates a small acoustic perturbation. 


ANALYSIS 

In the appendix the quasi one-dimensional equations governing acoustic propaga- 
tion through a variable-area flow duct are developed. The physical system is illus- 
trated in figure 1. An acoustic source located at x = 0 (j = 1) propagates sound 
upstream (in the positive x direction) against a steady, high subsonic Mach number, 
isentropic, compressible mean flow. The mean flow, travelling in the negative x 
direction, accelerates in the throat section as controlled by the area variation 
A(x) . 


Acoustics Equations 

The governing propagation equations, derived in the appendix, are as follows: 


u + u p + up] + -^(p u + u p + up) ^ = 0 


A^^s 


dx 


(1) 


3 


0 


(2) 


Su 

St 





P 




(3) 


where subscript s indicates a steady state mean flow quantity. All variables and 
coefficients in equations (1) to (3) have been nondimensionalized. (See the 
appendix.) Briefly, density has been nondimensionalized by p , velocity by c 
pressure by Pq^o' by multiplying by k, and time by multiplying by o), ^It 

should also be observed that the equations are nonlinear. 

Conservative form .- For the purpose of applying a finite difference scheme, 
equations (1 ) and (2) can be expressed in the following conservative form* 


St 


+ 



g 


0 


(4) 


Su 

St Sx 


0 


(5) 


where 


r(x,t) 


+ %P + up 


( 6 ) 


s(x, t) 


u 

u u + — 
s 2 


+ c 



(7) 


g(x,t) = -1/p u + u p + up\ ^ 
A^^s s^ dx 


( 8 ) 


Boundary conditions .- A harmonically varying source is assumed at x = 0 


p(0,t) = S sin t (t > 0) 


(9) 


where the amplitude S is calculated to give a specified source strength as dis- 
cussed later. A termination condition 


p(L,t) - u(L,t) = 0 


( 10 ) 
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is applied at x = L which, for a steady state, harmonically varying solution, 
reduces to the usual nonreflective impedance condition. Thus, the solutions for p 
and u are obtained from these equations, and p is substituted into equation (3) 
to get the pressure p. 


Numerical Method 

A fourth-order MacCormack differencing scheme (ref. 9) is used to generate a 
numerical solution of equations (4) to (8). The process has two steps applied at 
At/2 intervals. The first step involves a backward predictor and a forward cor- 
rector (BPFC). The second step uses a forward predictor with a backward corrector 
( FPBC ) . 

Consider equation (4) for the acoustic density and the discretization model 
illustrated in figure 1. Let p? be the n\imerical value of p(X‘,n At) f where 
j = 1, 2, ..., J and n = 0, 1, 2, ..., N. Assume that the values for have 

been calculated, and it is desired to calculate the corresponding values at the 
next At time step, that is, Pj 

BPFC phase .- In order to calculate p^"*"^ , the intermediate values are 

calculated for At/2 by using the backward predictor/forward corrector pRase of the 
differencing algorithm. The difference relations are (j = 1, 2, ..., J) 



^n _ At/2 _ 


Ax\^ 


8r 


n 


j“1 


3 - 2 ) 


At n 


( 11 ) 


pn+( 1/2) 




+ 


At/2 / (1) 

6 Ax\^ j 


8r 


( 1 ) 
j + 1 


+ 



At ^(1) 


( 12 ) 


The predictor quantities are denoted by a superscript (1), for example, Flux 

and force values based on predictor quantities are denoted in a similar manner, for 


example. 



and 



(See eqs. (6) and (8).) 


The predictor values Pj^^ at the next At/2 are computed from equation (11). 
Since j = 1/ 2, . . . , J, this requires flux values rQ and r_^ at the two extra 
grid points ahead of the source plane (j = 1). These are computed using third-order 
backward extrapolation as follows: 


n - n ^ n , 

r^ = 4r^ - 6r^ + 

0 1 2 


4r!J - 


(13) 


n 

1 




+ 4r!! - r"? 
2 3 


(14) 
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The predictor values ^ are used in equation (12) to calculate the corresponding 

corrector values n At + (At/2) time step. 

FPBC phase .- The next step is to use the corrector values p^+(1/2) 
n At + (At/2) in the forward predictor/backward corrector phase of the scheme to 
calculate p^ at n At + At* The required relations are (j = 1, 2, . . . , j) 




= + At/2 ^.^^n+(1/2) _ gj.n+d/2) ^ j.n+(1/2)'| n+(1/2) 


6 AxV j 


j+1 


■j+2 ) 




n+(1/2) 




(1) _ At/2 /'^^.(1) 


AxV^j 


- 

D-1 


+ r 


(1) 

j-2 


2 


(15) 


(16) 


The predictor values in this case require flux values rj^^ and ^j +2 the 

two extra grid points beyond the termination plane ( j = J) . These are calculated 
using third-order forward extrapolation of the variables p^'^^ "'/2) the BPFC 

step. These relations are as follows: ^ 


r 


( 1 ) 

J+1 



6r 


( 1 ) 

J-1 


+ 4r 


( 1 ) 

J-2 



(17) 


J+2 


= 4r 


( 1 ) 

J+1 


- 6r 


( 1 ) 


+ 4r 


( 1 ) 

J-1 


- r 


( 1 ) 

J-2 


(18) 


The same type of differencing relations (eqs. (11) to (18)) is applied to equa- 
tions (5) and (7) in order to compute the acoustic particle velocity u*?. Since it 
is a straightforward matter -to replace p^ by u^z these equations are^not given 

Once the solutions for p^ and u?? have been obtained, p^ is calculated 
from equation (3). ^ ^ 

Numerical boundary conditions .- In order to satisfy the source (eq. (9)) 
and termination (eq. (10)) boundary conditions, the characteristic variables 
p/Ps + u/Cg (for right-moving waves) and p/p - u/c (for left-moving waves) are 
used. For example, at any time t the instantaneous^density fluctuation at the 
source (x = 0) is given by equation (9). Thus, for the finite discretization, the 
source boundary condition is 


p^ = S sin t 
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The amplitude S is determined so that the pressure level at the source peaks at the 
desired decibel level (reference 2.0 x 10“^ N/m^). In the FPBC and BPFC sections of 
the difference algorithm, the quantity 


T 


1 




(19) 


is calculated at a previous time step where = p(0,t), u^ = u(0,t), p^^ = 

and Cg^ = Cg(0). Therefore, at any time t the source density is given by equa- 
tion (9), and the corresponding acoustic particle velocity is obtained from equa- 
tion (19); that is. 



At the termination section x = L, a nonref lective boundary condition is 
imposed. For a harmonic steady state this requires that the impedance of an outgoing 
wave be unity or, in other words, that p = u (eq. (10)). The required boundary 
conditions are obtained by calculating 


T 


2 



( 21 ) 


at a previous time step. This value for T 2 together with equation (10) and the 
linearized form of equation (3) yields the following boundary relations for p^ 
and Uj: 


P 


J 


T 

2 


c + 
sJ 




U 

j 


2 

c T 
sJ 2 


c + 
sJ 



) 


( 22 ) 


(23) 


where pj = p(L,t) / Uj = u(L,t), p^j = p^(L), and c^j = Cg(L). Note that the 
values of and T 2 are computed at time [n - (1/2)] At. The physical boundary 

conditions, however, are imposed at time n At. 
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DISCUSSION OF RESULTS 


The finite difference analysis discussed previously is applied in this section 
to a variable area duct containing high-speed subsonic flow and a high— intensity 
propagating acoustic wave. 


Duct Geometry 

The duct configuration is sketched in figure 1 and has an overall length of 
about 1 ,518 m« The two uniform sections on each end of the duct have a sguare cross 
section 0.0508 m on each side. The width of the duct is maintained constant 
(0.0508 m) as the height is gradually reduced to give a 2/1 area reduction at the 
throat (x/L = 0.50). Thus, the throat is rectangular (0.0508 x 0.0254 m) . 

The duct contour was designed so that a linear Mach number variation is achieved 
when the inlet (x/L = 0.25) and throat (x/L = 0.50) Mach numbers are 0.3 and 1,0, 
respectively. With this flow distribution, the one— dimensional gas dynamics egua— 
tions (eqs, (A1) to (A3)) can be solved to yield a closed form expression for the 
corresponding area variation A(x) . For a uniform section 


A(x) = 2.035 



(24) 


and for the variable area section 


A(x) 


1 r 

r + 1 

_ Y - 1 

G(x) 

2 

2 


G^(x) 


1 

Y-1 



(25) 


where G(x) is given by 


G(x) 


1 .34857 



(26) 


Sample Calculation 


In order to illustrate the finite difference algorithm, a calculation for the 
mean flow distribution shown in figure 2 is presented. Here the mean throat Mach 
number was arbitrarily set at 0.85. The distribution was then calculated by solving 
the one-dimensional gas dynamics equations using the area variation A(x) given by 
equations (24) to (26), The Mach number increases from about 0,28 at x/L = 0.25 
to 0.85 at the throat and then drops to about 0.28 at x/L = 0.75. Hence, the cal- 
culated mean flow is symmetric about the throat (x/L = 0.50). For the acoustic data 
which follow, a negative throat Mach number is indicated (M = -0.85) to emphasize 
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that the flow is in the negative x direction, opposite to the propagation of acous- 
tic waves from the source. This corresponds to the aeroengine inlet situation. 

Figure 3 shows the acoustic pressure distribution calculated by the nonlinear 
finite difference algorithm (eqs. (11) to (23)) for a 1 35-dB acoustic source oscil- 
lating at 500 Hz. For comparison, the results of a linear finite element calculation 
are also shown. The finite element algorithm is based on the Galerkin/weighted- 
residual method using one-dimensional, isoparametric, Hermite elements of third order 
(see ref. 11, for example). 

Owing to the nonlinear character of the finite difference model, the solid curve 
in figure 3 is actually the fundamental (500 Hz) component (sin t) of the acoustic 
pressure. The two calculations, nonlinear finite difference and linear finite ele- 
ment, are in close agreement. This is unexpectedly true in the high-intensity throat 
region. Spectral data are also presented in table I for the fundamental and first 
three harmonics of the nonlinear finite difference pressure calculation. Comparison 
of these data indicates that a considerable amount of energy is shifted into the 
higher harmonics. For example, at x/L = 1.00, the amount of energy contained by the 
first harmonic (sin 2t) is about 16.5 percent of that contained in the fundamental. 

A plot of the acoustic particle velocity is shown in figure 4, and the acoustic 
streaming velocity component is presented in figure 5. With a linear acoustics model 
driven sinusoidally, the acoustic velocity would have zero mean. With nonlinearities 
present, however, the acoustic velocity exhibits, in addition to higher harmonics, a 
nonzero mean component termed the streaming velocity. A nonzero streaming velocity 
indicates a transfer of acoustical energy to the mean flow. As shown in figure 5, 
the acoustic streaming velocity for the present case peaks in the throat region, but 
compared with the mean flow (fig. 2) , is quite small. 

Figure 6 shows time histories of the acoustic pressure at five axial locations 
after steady state was achieved by the finite difference algorithm. Figure 6(a), 
shown for reference, is the source pressure time history, which varies sinusoidally 
at 500 Hz (sin t) . The peak amplitude corresponds to a level of 135 dB. At station 
x/L = 0.25 (fig. 6(b)), which is located at the beginning of the area variation, the 
pressure response is still very close to a pure sine wave. For example, table I 
shows a first harmonic level of 107.1 dB or about 27.7 dB lower than the fundamental 
(134.8 dB). 

At the throat (fig. 6(c), x/L = 0.50), the pressure shows appreciable distor- 
tion from a pure sine wave. Here the pressure peaks at about 150.1 dB (table I and 
fig. 3). The results of the linear finite element indicate a peak of approximately 
150.0 dB. The first harmonic (137.4 dB) in this instance is only about 12.7 dB lower 
than the fundamental (150.1 dB). It thus appears, by virtue of the increased har- 
monic content, that the nonlinearities are becoming more important, although the 
linear solution still adequately predicts the peak level in the throat region. 

Further upstream at the end of the area variation (x/L = 0.75), as shown in 
figure 6(d), considerable harmonic distortion remains although the peak level has 
dropped rapidly to about 134.6 dB. Here the first harmonic (126.5 dB) is only about 
8.1 dB lower than the fundamental (134.6 dB). Finally, at the termination 
(fig. 6(e), x/L = 1.00), about the same picture prevails, although the beginning of 
a steepening in the wave is noticeable. 
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This wave steepening behavior is examined further in figure 7, where the source 
strength is varied from 134 to 137 dB. Here a very definite nonlinear steepening 
appears in the pressure time history at the termination. This indicates that by 
increasing the source strength, the acoustic response approaches a shock condition. 


Comparisons With Asymptotic Theory 

As discussed previously, the finite difference algorithm was compared with a 
linear finite element solution for the classical small amplitude linear acoustics 
case. Here a check is made for a typical nonlinear situation by comparing the finite 
difference solution with the asymptotic theory of Callegari and Myers (ref. 10). The 
typical nonlinear situation arises at higher sound pressure levels for throat Mach 
numbers approaching unity. In this theory (l - |m^|) is taken as the (small) per- 
turbation parameter. 

A slightly different geometry was used for these calculations. As shown in 
figure 8, the constant area entrance and termination sections were eliminated, and 
the area variation is given by equation (25) where 


1 

2 


G(x) = 1 - 0.91287 — 

L 


(o < f < l) 


(27) 


For this geometry, the area reduction ratio is about 1.32/1, so that for an entrance 
Mach number of about 0.51, the flow accelerates linearly to a sonic velocity at the 
throat. Hence, this area variation provides a much more gradual choking of the flow 
than does the configuration (eqs. (25) and (26)) used in the earlier calculations. 

Results are presented in figures 9 to 1 1 for three combinations of source 
strength and throat Mach number. In each case, the acoustic velocity divided by 
(’’ “ l^tl ^ “ 0.75L and x = L is shown after steady state had been achieved 

by the differencing algorithm. For the cases presented, the quantity (l - |m^|) 2 
is of the order of 10 or less. Hence, the calculations are accurate to at least 

three decimal places. The solid lines in figures 9 to 1 1 indicate the finite differ- 
ence solution. The circular and square symbols indicate the calculated values based 
on the asymptotic theory of Callegari and Myers (ref. 10). 

The comparisons for a throat Mach number of M^ = -0.90 are graphed in fig- 
ures 9 and 10. The strengths, for an equivalent pressure source located at x = 0 
(fig. 8), are 130 and 140 dB, respectively. In both cases, no shocks were predicted 
or, for that matter, were expected, and the finite difference calculation agrees very 
well with the asymptotic theory. At 140 dB, as shown in figure 10, a severely dis- 
torted time history occurs at x/L = 0.75, and the agreement is excellent. 

Figure 11 illustrates a similar comparison for a 130-dB source when the throat 
Mach number is increased to M^ = -0.95. This situation is very close to shock for- 
mation (ref. 10). Here the comparison is still very good. However, there are some 
minor differences between the asymptotic theory and the finite difference algorithm 
near t « ti /2 for x/L = 0.75. 
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For this latter case (fig. 11), a further increase in the throat Mach number (or 
an increase in the source strength) will precipitate a shock and a resulting loss of 
acoustical energy. For weak shocks, the energy relation given by equation (3) might 
suffice. Stronger shocks, manifested by appreciable loss in acoustical energy and 
increase in entropy, will require that the general energy equation (ref. 12) be sub- 
stituted in lieu of equation (3). In both instances a "shock" jump condition must be 
added to the finite difference code. 

Finally, it should be noted that the asymptotic theory of reference 10 is 
extremely accurate when (l - |m^|) approaches zero; that is, when the throat Mach 
number approaches unity. It is precisely these situations which present the most 
computational difficulty for the finite difference algorithm, since axial wavelengths 
become increasingly small, which necessitates extremely fine time/space meshes. 
Therefore, the comparisons presented in figures 9 and 10 are thought to represent a 
stringent test of the differencing procedure to simulate the nonlinear acoustics 
case. The finite difference algorithm is more effective in handling the calculations 
for larger values of (l - With these situations the asymptotic theory is less 

accurate. Hence, these two independent methods nicely complement one another. 


CONCLUDING REMARKS 

An explicit finite difference algorithm based on a MacCormack scheme was used to 
solve numerically the one-dimensional nonlinear acoustic equations for a variable 
area duct. The procedure iterated over the time/space discretization grid until a 
steady state solution was achieved. Sample calculations were presented for several 
combinations of area reduction ratio, throat Mach number, and source sound pressure 
level. 

The principal effect of the nonlinearities was to shift acoustical energy to 
higher harmonics. With increased source strengths, wave steepening was observed. 

This phenomenon suggests that the acoustical response may approach a shock behavior 
at a higher sound pressure level as the throat Mach number approaches sonic value. 
Where applicable, comparisons were made with the calculations from a linear finite 
element algorithm. On a peak level basis, good agreement between the nonlinear 
finite difference and linear finite element solutions was observed, even though a 
peak sound pressure level of about 150 dB occurred in the throat region. Nonlinear, 
steady state waveform solutions were shown to be in excellent agreement with a non- 
linear asymptotic theory. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton , VA 23665 
June 9, 1983 



APPENDIX 


DERIVATION OF NONLINEAR ACOUSTIC EQUATIONS 

In the following, the one-dimensional compressible fluid continuum equations are 
used to derive the governing acoustic propagation relations. 


Fluid Continuum Equations 

The one-dimensional compressible flow equations for an isentropic fluid may be 
expressed in terms of the density p and velocity u (ref. 12). The equation of 
continuity is 


_ A# 

•^ + 4 ( 5 =) 0 

9t 9x A dx 


(A1 ) 


and the balance of momentum is given by 


^ 9 /u^ 

5t 9x\ ^ 




(A2) 


The pressure is determined from 
^ Po 

P = — P (A3) 

^o 

In the above equations p^ and are the ambient pressure and density, 

respective ly . 

Acoustic Perturbations 
In equations (A1) to (A3) let 

p(x,t) = pg(x) + e p(x,t) 

u(x,t) = Ug(x) + e u(x,t) 

p(x,t) = Pg(x) + e p(x,t) 

1 2 



APPENDIX 


where a caret ('") over a symbol indicates a small acoustic perturbation about the 
basic steady flow state, denoted by a subscript s. The steady state flow equations 
are 



+ Ps^s dA 
A dx 


0 


(A5) 


5x 



+ 




(A6) 


~ Po 

Pc? = P 

^S y f'Q 

P 

^O 


(A7) 


Substituting equations (A4) into equations (A1) to (A3) and expanding and then sub- 
tracting steady state terms (eqs. (A5) to (A7)) gives the final nonlinear acoustic 
equations 


9t 



M A 

UgP + 


c55) + 4(i 


A A \ 

Ps" + ^sP + e^p) 


dA 

dx 


0 


(A8) 


du 

at 





(A9) 


The acoustic pressure is given by the relation 


P = 



(A10) 


where the local speed of sound is given by 


~2 

^s 


TfPs 

Ps 


YPo 




(All) 
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APPENDIX 


The parameter e denotes a nonlinear acoustic term. With e = 0, equations (A8) 
to (A1 0) give the usual one-dimensional linear acoustic solution. Whereas with 
e = 1, a nonlinear solution is obtained. 


Nondimensional Acoustic Equations 

Equations (A8) to (A10) are nondimensionalized by using the following 
transformations : 


A ^ 


p{x,t) = 

Po 


u(x,t) = 


^ 

u(x^ t) 


p(x, t) 


p(x, t) 
2 

P ^ 
'^o o 


X = kx 


k 



Ug(x) 

Ug(x) = 

<=0 


> 


A(x) 


A(x) 

At 


pg(x) 


Po 


t = (jL)t 


C 


s 




(A1 2) 


14 



APPENDIX 


Substituting equations (A12) into equations (A8) to (A10) yields the nonlinear acous- 
tic equations in a final nondimensional form 


5t 


+ 


6x 



u p 


e»p) <-x(ps“ 


+ u D 
S ^ 


+ eup^ 


dA 

dx 


0 


(A1 3) 


5u 

at 



+ 



0 


(A14) 


P 



+ e 


I— I 



(A1 5) 


2 

c 

s 



(A1 6) 


In terms of these nondimensional variables, the instantaneous local Mach number M 
is 


M 


+ eu 


izl 

(ps + ep) ^ 


(A17) 
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TABLE I.- ACOUSTIC PRESSURE FOR 135-dB SOURCE 
[M^ =0.85; f = 500 Hz] 


X 

Nonlinear finite 

difference, dB 

Linear 

finite 

L 

Fundamental 

1st harmonic 

2d harmonic 

3d harmonic 

element 

dB 

0.00 

.25 

135.0 

134.8 


84.7 

90.3 

1 35.0 
134.7 

.50 

150.1 


127.9 

102.2 

150.0 

.75 

134.6 


116.9 

105.6 

135.2 

1 .00 

134.7 

126.9 

117.1 

105.4 

135.3 














Velocity, u Sound pressure level, dB 


O 


Nonlinear finite difference 
Linear finite element 



Figure 3.- Acoustic pressure. 135-dB source; = -0.85; f = 500 Hz. 



Figure 4.- Acoustic particle velocity. 


135-dB source; = -0.85; f 


500 Hz. 
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